HEXAGONAL LATTICE POINTS ON CIRCLES 



OSCAR MARMON 

Abstract. We study the angular distribution of points in the 
hexagonal lattice (i.e., ^[ i+v^ ]) lying on a circle centered at the 
origin. We prove that the angles are equidistributed on average, 
and show that the discrepancy is quite small for almost all circles. 
Equidistribution on average is expressed in terms of cancellation 
in exponential sums. We introduce Hecke L-functions and investi- 
gate their analytic properties in order to derive estimates on sums 
of Hecke characters. Using a version of the Halberstam-Richert 
inequality, these estimates then yield the desired results for the 
exponential sums. 

An interesting consequence of these bounds is that the discrete 
velocity model (DVM) for the Boltzmann equation is consistent 
when using a hexagonal lattice. 
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1. Introduction 

1.1. The Hexagonal Lattice. We will study the properties of a lat- 
tice in that is spanned by the vectors (1, 0) and Q, see Figure 
d We are interested in the following two questions: given a circle with 
radius r, centered at the origin, whose perimeter contains at least one 
lattice point, 

(1) how many lattice points lie on the circle and 

(2) how are these distributed around the circle? 

Why are these questions interesting? Well, aside from the intrinsic 
number theoretic and geometric importance, there is, perhaps surpris- 
ingly, an application to the Boltzmann equation in the kinetic theory 
of gases. 

1.2. The Boltzmann Equation. Under certain hypotheses, the be- 
havior of a gas is described by the Boltzmann Equation: 



(1) 



dt dx 



QifJ)- 



Here f{x,v,t) : x M'^ x M"*" is the phase space density of the 

gas, that is, it describes the expected mass density of the gas at time 




-3 -2 -1 



Figure 1. The hexagonal lattice 
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25i— r 




-30 -20 -10 10 20 30 



Figure 2. Lattice points on the circle with radius r = 21 



t in the point {x, v) in six-dimensional phase-space, where a; G M'^ is a 
position in space and f G is a velocity. 

If we replace 3 with any spatial dimension d, the Boltzmann equation, 
although losing its physical sense, is still of mathematical interest. We 
will consider the case d = 2. Q{f, f) is called the collision term and is 
given in the 2-dimensional case by 

(2) Q{f,f){v) = J (j{f{v')f{v:)-f{v)f{v,))q{\w\,cose)de\dv,, 

where are the velocities after the collision and v',vl those before. 
Moreover, cos^) is a quantity describing the probability that two 

particles with relative velocity 2w, i.e. 

V — V^: 



collide with relative deflection angle 6. We have the relations 




v' = i(f + f^,) + \w\u, 
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where m is a unit vector with u = {cos 6, sin 6), which follow from the 
laws of conservation of energy and momentum. In other words, v' and 
are endpoints of a diameter of a circle with radius \w\, centered at 
{v + v,)/2. 

In a discrete velocity model (DVM), one considers a discrete set of 
possible velocities. (See PP for more details.) Here we will consider a 
DVM where the velocities belong to the set hh, where h > and L is 
some lattice in M?. Set 



5eL 



so that f'^ f as h 0, in some suitable sense. We want to prove 
the consistency of the DVM, i.e. the property that 



Q{f, f){v) for all v e hh, as h ^ 0. 



g,iw,e) = (fiv')fiv:) - fiv)fiv,))qi\wlcose). 



(4) 
Let 



Then, by a change of variables 
(5) Q{f,f){v)- 




gv{w, 9)d9 dw. 



When we discretize the integral in (0), the outer integral turns into a 
sum over all lattice points w = h(, ( & h, while for the inner integral we 
get a sum over all u such that v' and v'^ belong to hh. A necessary and 
sufficient condition for this is that m = |||, where ^ G L and |^| = 
Thus we get 



1 



9v{hC,aTg^) 



V 



ICNICI 



where is the area of a fundamental region of L and 

rdn) := ^ e L; \^\' = n}. 

It is clear that gv{w,6) is a 27r-periodic function of 6, and so can 
be expressed as a Fourier series. It turns out that, assuming certain 
regularity conditions on gy{w, 9), a sufficient condition for (jlj) to hold is 
that we have sufficient nontrivial cancellation in the exponential sums 

S{n,A):= e^^'^-'s^'^) 

AtSL 
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as n grows, ii A & Z, A ^ 0. 

Fainsilber, Kurlberg and Wennberg proved ^ that if we choose L to be 
Z^, then (jU holds provided gy{w,9) is a C^-function. This is done by 
identifying with the ring of Gaussian integers, Z[z] C C, and apply- 
ing number theoretic methods. The results of this paper indicate that 
we can equally well consider the hexagonal lattice defined above, corre- 
sponding to the ring of integers in the algebraic number field Q(\/— 3). 
Similar methods to the ones used here should give the same results for 
every imaginary quadratic number field Q{^/—d) with class number 1 
(and thus unique factorization into irreducibles) . There are exactly 9 
such fields: Q(\/— d) has unique factorization if and only if d takes one 
of the values 1, 2, 3, 7, 11, 19, 43, 67 or 163 (see [Hj, Ch.l3,§l). 
Moreover, it is reasonable to expect (jlj) to hold for any lattice L such 
that 

limsuprL(n) = oo. 

n— >oo 

For dimensions li > 3 (i.e., including the physically relevant case d = 
3), the consistency of a DVM based on the lattice Z'^ was proven by 
Bobylev, Palczewski and Schneider [T], using deep number theoretic 
results. 

1.3. A Number Theoretic Point of View. To answer the questions 
on Page El we will view the lattice not as a subset of but of C, and 
use number theoretic methods. We identify the lattice with the set 

Zfcj] = {a + buj ] a,b & Z}, 

where 



This is in fact the ring of integers D in the algebraic number field 
K := Q(a/— 3). Since K is a. quadratic number field, the norm N{a) 
of an element a & K is just the squared complex modulus Hence 
question ((H) above can be rephrased as: 

How many elements of O are there with a given norm? 

Let a, 6 G Z, so that a + buj E O. Then we have 

N{a + buj) = (a + buj){a + bZJ) = + (cu + uJ)ab + touJb'^ = + ab + 6^, 

so that yet another formulation of the first question would be: 

How many integer solutions (a, b) exist to the equation 
a? + ab + b^ = R for a given R e Z7 
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This is a classical question, and the answer we present is by no means 
new (see for example , Ch. 12.9). The norm of a quadratic number 
field induces a quadratic form, in our case 

Q{x,y) =x^ + xy + y'^. 

The quantity we seek is then 

rQ{n) := #{(a, 6) G g(a, 6) = n}, 

the number of representations of n by the form Q. The key point is, 
that each lattice point a+hu on the circle with radius r corresponds to a 
representaion of as a product of two elements a + hu and a + 6uJ in D. 
Fortunately, D is a principal ideal domain, hence a unique factorization 
domain, so all such representations can easily be found by factorizing 

into rational primes. 
For every rational prime p, one of the following three cases is true: 

• p = TT -W where vr is a prime in D. We say that p is split in D. 

• p = n"^ where tt is a prime in O. We say that p is ramified in 
D. 

• p remains prime in £). We say that p is inert in £). 

These factorizations are unique up to multiplication by a unit. The 
unit group of Z[lj] is 

U = {±1, ±LU, ±LU^} = {UJ, W^ . . . , UJ^}. 

Moreover, all primes in 1j[uj] can be found in this way. For odd p it 
is known that (see for example jT7j, Ch. 13, §1) p is split if (^^j = 1, 

ramified if (^^^ = 0, and inert if (^^^ = —1. Moreover it is easily 
seen that 2 is inert, since the equation + a6 + 6^ = 2 has no integer 
solutions. We thus have: 

Proposition 1. Let p be a rational prime. Then 

(1) p is split ijf p = 1 (mod 3). 

(2) p is inert iff p = 2 (mod 3). 

(3) p is ramified iff p = 3. 

We see that to each rational prime p with p = 1 (mod 3) corresponds 
a unique prime Hp = ^e*^p in Z[uj], if we choose 9p to lie in the interval 
[0, |]. Indeed, for every tt, one of the numbers un, . . . , u^n, uW, . . . , u^W 
must lie in this interval. For p=3, analogously, if we set 
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we have 3 = nsW-^. We note however that ¥3 = uJtts, so that 3 = cUttI, 
confirming that 3 is ramified. 

Now let n G Z and suppose that n has the rational prime factorization 

Pi ■■■Pk Qi 

where pi, . . . ,pk = 1 (mod 3) and qi, . . . ,qi = 2 (mod 3). For each pi 
we choose the unique prime TTj = y/pe^^p^. Then the complete factor- 
ization of n in Z[co'] is 

for some a G N. It follows that if we have the factorization 

(6) n= {A + Buj){A + BU) = + AB + B\ 

then the multiplicity of vTj in A + Buj can be any integer m with < 
m < ai, and the multiplicity of vfj must then be — m. 
A necessary and sufficient condition for n to be representable in the 
form (jH)) is that all the (3i are even numbers. However, the qi have no 
effect on the number of solutions nor on the arguments of the solutions 
A + Blj. 

The power a of 3 in n, increases the argument of each solution by a 
multiple of ^, but has no effect on the number of solutions, since the 
prime factors of 3 are associated. Finally, A + Buj can be multiplied 
by any unit. This yields that the number of solutions to (jHl) is 

(7) rQ{n) = 6- J] + 

P^ = l (3) 

Example. The circle in Figure |21 has radius r = 21. The lattice points 
on the circle therefore correspond to factorizations of 

in Z[ct;]. Here 

3 = {iV3){-iV3), 7= (3-cj)(3-cJ). 

If /i = A + Bu lies on the circle, so that 21^ = fj,Ji, then /x is one of the 
three numbers 

3(3-cj)2, 21 or 3{3-ijf 
multiplied by any of the six units, giving 18 lattice points. 

We have answered the first question on page El and consequently 
devote the rest of the paper to the second one. 
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1.4. Equidistribution through Exponential Sums. We want to 
study the distribution of hexagonal lattice points on circles. What we 
would like to conclude is that the arguments of all the lattice points of 
a given circle is in some sense evenly distibuted in the interval [0, 2it). 
Let us define equidistribution of a sequence in the following intuitive 
way: 

Definition. A sequence {xn}nei+ ^ [^j B) is said to be equidistrihuted 
if 

^.^ #{n <N]Xne {a,h)} ^ b-a 
N^oo N ~ B - A 

for all A < a < 6 < B. 

Equidistribution can be formulated in terms of cancellation in expo- 
nential sums as follows: 

Theorem 2 (Weyl's Criterion). The sequence {xn} is equidistrihuted 
in [0, 27i) if and only if 

1 ^ 

hm — y e-*^^" = 

n=l 

for all integers k ^ 0. 

Proof. The proof is based upon the observation that a characteristic 
function of an interval can be approximated by trigonometric polyno- 
mials. (See m, Ch. 1, Th. 2.1.) □ 

To what extent are the arguments of hexagonal lattice points on 
circles equidistributed? In Section |S] we will show that there exist arbi- 
trarily large circles with arbitrarily poorly distributed lattice points, so 
we do not have equidistribution for all circles. We can however prove 
that lattice points on circles are equidistributed on average. Inspired 
by Weyl's Criterion, we introduce exponential sums: 

Definition. Let n G N, A G Z \ {0}. Define 

Sin, A) = e'^-"'^". 

|/.|2=n 

Then the main result of this paper, states that the arguments of 
hexagonal lattice points on circles are equidistributed on average in 
the following sense: 
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Theorem 3. If 6 \ A then S{n,A) =0. If A ^ Q and Q \ A we have 
for every 5>1— as x^oo and A = 0<5(eV^^°^), 

-J2\S{n,A)\<^s (logx)-^ 

n<x 

One can interpret this resuh in the foUowing way: the trivial estimate 
for En<x \Sin,A)\ is 

j:\S{n,A)\<J2rQ{n)^nx. 

n<x n<x 

By Theorem ini however, the sum is o{x), ergo the terms are in some 
sense smaUer than expected. This is enough for the apphcation to the 
Bohzmann Equation. 

We will also show that the arguments of the hexagonal primes, or 
equivalently, prime ideals, are equidistributed. We can define a unique 
argument of every prime ideal in the following way: 

Definition. Let p be a prime ideal. If p = (a) there is a unique a' 
among the associates of a satisfying — f < arga' < |. We then define 
9p = arga'. 

We will then prove 

Theorem 4. The sequence {Op}, where p ranges over all prime ideals 
of O, ordered by growing A^(p), is equidistributed in [— |, |). 

1.5. Outline of what to come. In Section |21 we will introduce two 
number theoretic concepts related to the exponential sums described 
above - Hecke characters and Hecke L-functions. The definition of 
the Hecke characters is in general much more complicated than in this 
special case, where matters are simplified by the unique factorization in 
the ring of integers of K. The main result of Section |21 is the functional 
equation for the Hecke L-functions, a special case of the ones derived 
by Hecke for general algebraic number fields in [H]. 

In Section El we will use analytic methods to estimate a sum of Hecke 
characters over prime ideals. The methods are essentially those that 
are used in the classical proofs of the Prime Number Theorem and 
the Prime Number Theorem for Arithmetic Progressions (with error 
terms). In his paper 0, Kubilius proved these results for the Gaussian 
number field Q{i). 

In SectionEJwe will prove equidistribution on average of lattice points 
on circles, using a lemma on mean values of multiplicative functions. 
The method is similar to the one used in the papers by Katai-Kornyei 
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pnj and Erdos-Hall [121 concerning the distribution of Gaussian inte- 
gers. 

Section |31 contains a theorem about the density of prime ideals in 
circle sectors, and a construction of circles with many points that are 
poorly distributed. 

Finally, in Section IHl we give an estimate for the discrepancy, which 
is another measure of equidistribution. 

2. Characters and L-functions 

2.1. Definitions and Basic Properties. Let K = Q(v^— 3) and let 
O = 1j[uj], the ring of integers in K. If a is an ideal in D, we will write 

a<0. 

A Hecke character (modulo (1)) of the number field i^' is a function 
^6a fj-QYa the ideals of O to the unit circle defined by: 

This is well defined since every ideal in O is principal, and since (/i) = 

(z/) implies ^ = eu, where = 1. 

The Hecke L-function L{s, x^") is defined by 

^y^^X ) A^/ I |2s+6a c |,,|2s+6a' SXHyb)^!. 

(Following the conventional notation, we use the complex variable s = 
a + it.) 

Remark. When a = 0, the L-function reduces to the Dedekind zeta 
function 

We must of course prove that the series defining the L-functions 
converge. First of all we note that the series 

1 

where p runs through all prime ideals in D, converges for a > 1. Indeed: 



E 



EtttV < 2y]p"" < 2 Vn-" < oo. 



Now we prove 
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Proposition 5. The series 

sr- X'°(a) 

converges absolutely for a > 1 and uniformly for a > 1 + 6 > 1, making 
L(s, x^") analytic function for a > 1. Furthermore, for a > 1, we 
have the following Euler Product representation: 



p<0 



N{py 



where the infinite product is absolutely convergent. 
Proof. The general factor of the right side of (jH} is 



This is bounded by 

since A^(p) > 2 for all prime ideals p. Thus, by the above remark the 
product in © is absolutely convergent for a > 1. 
Now, 

V Nipy J [ ^ NipY + iv(p)- + • • 

(9) ^E^^"^^^ 



N{a)<x Af(a)> 

where the star indicates that a runs through those ideals whose divisors 
all have norm < x. In particular, letting a = and s = a > 1 in (jH)), 
we see that 

-1 



N{a) 

establishing the convergence of 
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Also, by © 



< E 



when X ^ oo, which completes the proof. □ 
Proposition 6. For Re(s) > 1 

Proof. By (jH)) we have, for cr > 1 

-logL(s,x^ ) = 2^ log 1 - — — ]=-2^2^ 



Differentiation yields 



L{s,x'^) V^i ^(P) 



Remark. From now on we assume that a 7^ 0. 



□ 



2.2. A Theta Formula. The analytic continuation of the L-function 
to the whole plane will be proved through a transformation formula 
for a so called theta-function. We will give the theta formula that was 
proved by Hecke fSj, in the particular case of the field K = Q(-\/— 3). 
The idea behind the proof is to use two-dimensional Poisson Summa- 
tion. We shall recall some Fourier theory. For variables in it is 
always to be understood that x = {xi,X2), y = (1/1,1/2) and so on. 

Definition. A function / G C°°(]R^) is called a Schwartz function if it 
approaches zero faster than any inverse power of CC clS |x| ^ 00 5 clS do 
all of its derivatives. 

Definition. We let (■, ■) be the standard scalar product in M^: 

{x,y) := xiyi + X2y2- 

Definition. For a Schwartz function / we define the Fourier transform 
/by 

f{y) = I f{x)e'-'^^'yUx. 
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Theorem 7 (Poisson Summation Formula) . Let f be a Schwartz func- 
tion. Then for every x we have 

Proof See Lang ({H|, XIII,§1) □ 
The Fourier transform of the function 

is particularly simple - it is easily seen that we have 
(10) h = h. 

Definition. Let / be a Schwartz function, B a non-singular real ma- 
trix. Define 

fnix) = f{Bx). 

Obviously, fs is then also a Schwartz function. The Fourier trans- 
form is easily found: 

Lemma 8. We have 

My) = j^^f{{B-Yy), 

where \B\ is the absolute value of the determinant of B. 
Proof. We have 

fB{y) = I f{Bx)e-^^''^^^y^dx. 
By the change of variables z = Bx we get 

My) = ^ / fiz)e-'-^^''-'^^y^dz 

= 1 f f{z)e-'-^(^'^^-'^^yUz 
\B\ J 

= j^^fiiB^Yy). 

□ 

Definition. Let Q{x) be a positive definite quadratic form given by a 
matrix A (that is, Q{x) = {Ax,x)). We define the Schwartz function 
fqix) by 

fgix) = e-'^Q^^). 
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Lemma 9. We have 



where Q'{x) = {A ^x,x). 

Proof. Since A is positive definite, there exists a real matrix B such 
that = A. Then Q{x) = {Ax, x) = {Bx, Bx), so that we have 

fgix) = heix), 

where h(x) = e^^^^'^K Thus, by Lemma |H1 and (|lUj). 

Uy) = hsiy) = T^h{B-'y) = -^h{B-'y) = -^e--Q'^y\ 

II vl^l vl^l 

□ 

Now we prove a formula from which the desired theta formula will 
follow. 

Theorem 10. Let Xi,X2 be real variables and define 

^^^^ {ui=Xi+UJX2, 

]U2 = Xi+ ZJX2. 

Moreover, lett be strictly positive. Then we have the following formula: 

(12) 6"^''*^''+"^)'^''+""^ = —5— e"^''''''^^^'^"''"'''"^ 

Remark. For real xi,X2 we have ui = U2- Later we will extend this 
result to complex xi,X2 in which case this relation does not hold in 
general. 

Proof. Define the positive quadratic form Q{y) by 

Q{y) = 2t\yi + 7/2^1' = 2t{yl + yiy2 + yl). 
Then Q{y) = {Ay,y), where 

Since 
we get 

Q'{y) = {A-^y, y) = ^(y^ - y^y^ + y^). 
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Now let US define 
Then we have 

= e-2"*l"i+"2^+^i+^2^l' = ^ e-"^(^+'")= ^ fQ{x + m). 

Using Poisson Summation and Lemma El we get 

F{x)= J2 /QMe^'^*^™'^^ = J2 e-"^'("^)+2"^<™'^\ 

Solving (fTTj) for Xi,X2 yields 



xi = ■:y^{ujui - UJU2] 

X2 = 75 (-^-^1 +^^2)- 



Thus we get 



(m, x) = rriiXi + 7712X2 = — ^(^Mi(miCt; — 7712) — U2{miuj — m2)j. 

Put u = miUJ — m2. Then if m runs through Z^, runs through D. 
Moreover 



Q'{m) = ^("ii - mim2 + ml) = 



and 



(m, x) = —r=[h>ui — 77^2). 
v3 



We conclude that 



as stated. □ 

We will now extend this result to complex numbers xi, X2. As before 
we let 

Ml = Xi + ^01X2, 
U2 = Xi+ LJX2. 

Note that when we allow arbitrary complex values for xi and X2, Ui 
and U2 become independent complex variables. 

Proposition 11. The series 

F{x) = J2 e~2-t(M+-i)(M+-2) and G'(x) = e-^l"l'+^(""^-""^) 
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are absolutely convergent for all {xi,X2) G C^. Moreover, they converge 
uniformly in the region Q^i = {{xi,X2) G : max(|a;i|, \x2\) < R} 
for every R> 0. 

Proof. Let 

Vf,{x) = 27lt{fl + + U2), 

so that 

F{x) = ^ e-^"("). 

Expanding 

V^{x) = 27rt(|/i|^ + yUM2 + 'pui + U1U2), 
we see that, in VLr, 

Re(P^(x)) > 27rt(|/i|2 - 4|/i|i? - AR^) > 
Thus, for all but a finite number of /i, 

and e~^''^'^ clearly converges, so the uniform convergence of the 
series F{x) follows. Moreover, 

G{x) = J2 e~®^^"\ 

where 

27r 2 27r 
" 3t - --/|(z/U2 - z^ni). 

In fi/j we have 

Re(Q,(x)) > — |z/p - -^|z/| > |z/|2, 

so the uniform convergence of G{x) follows by an analogous argument. 

□ 

Thus both sides of fll2|l define entire functions of the complex vari- 
ables Xi,X2- Since they agree for real Xi,X2, they must be equal. Thus 
we have proven 

Theorem 12. Fort > and arbitrary complex numbers Ui,U2 we have 

(13) g-27rt(At+ni)(77+«2) ^ _J_ ^ fl {i^«2-i^«i) _ 

From Theorem El we derive our theta formula: 
Definition. Let 
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Theorem 13. We have 

e{t,a) = t-^-^''e{j,a) 

Proof. We let p G D be arbitrary. Furthermore we introduce the vari- 
able z, and set 

Ui = p U2 = Z + p 

Now (fT5|) is equivalent to 



(14) ^expj -27rt(|/i + p|2 + 2;(^ + p))| 

1 V- f 27r, ,2 27r,_, _, ,1 

We differentiate p4|) 6a times with respect to z. This yields 

(15) ^(/i + p)6"exp| -27rt|/i + p|2 + ;z(/i + p)| 



^ ^ ^I76"exp|-^|z/p + ^(l7(2; + p)-z/p)|. 



Setting 2 = 0, we get 

(16) ^(/i + p)^"exp| -27rt|/i + prt 

Put t = As yU ranges over O, so does p + p, since p G O. Thus the 
left side of (fTB|) equals 

Y^p'^e-fs^^'^" = 9{T,a). 

We note that 

1 2 

-^(iTp - z/p) = ~i—lm{up) 

But Im(z/p) = for some m G Z, and thus we can neglect the term 
-^(/Ip — pp) in the exponent of the right side of (|T!)|l . Therefore the 
right side equals 

since when u runs through D, so does I^. This finishes the proof. □ 
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2.3. The Functional Equation. We will now see how the L-functions 
can be extended to the whole complex plane. 

Definition. Let 

Theorem 14. is entire and satisfies the functional equation 

Proof. Since L(s,x^'*) = L{s,x^'^~°^) we can assume that a is positive. 
For fi E K and a > 1 we have 

oo 


and hence 



r(s + 3a)^^= //i^'^e-^l'^l't^+^a^irft 



\V3 
Thus 



r, \ s+3a oo 



(/Tj^ s+3a oo 

oo 



(17) 

~ 6 



For Re(s) large enough, (fTTjl is absolutely convergent, and hence we 
can change the order of integration and summation to get 

(r^\ s+3a oo 

(18) 



6 



The right side of (fTHj) is 



-.1 ^ oo 

6 J 6 J 
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Using Theorem El this is 



1 




6 



1 




1 1 
Thus we have proven (for Re(s) large enough) 




1 



Iv. 



But since this integral converges absolutely for all s, ()19|) represents an 
analytic continuation of C,{s, x^") to the whole plane. The functional 
equation also follows, since the right side of (jl9|) remains unchanged 



From this theorem we see that L{s, x^"") is an entire function. Using 
our knowledge of the Gamma function, we also deduce that L{s, x^"") 
must have zeros at s = — 3|a|, — 1 — 3|a|, — 2 — 3|a|, .... These are called 
the trivial zeros. Moreover, for cr > 1 it can be seen from the Euler 
product that L{s, x^") ^ 0. 

There is however an infinite number of non-trivial zeros in the so called 
critical strip, < a < 1. It is generally believed that all of them lie on 
the line Re(s) = |. This statement is part of the generalized Riemann 
Hypothesis, and if we assume it, the estimates in SectionElcan be made 
much sharper. In the next section, however, we will unconditionally 
narrow down the region in which the non-trivial zeros can appear. 

Corollary 15. In the strip — ^ < cr < 4 



when we replace s by 1 — s. 



□ 



L{s,xn = he''' 
for some positive constants ki, k2 (depending on a). 
Proof. By (IT^ . 




1 



and hence 
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Stirling's Formula states that, in the angular region — tt + 6 < args < 
TT + 6 for any fixed 5 > 0, we have as |s| ^ cxd 



logr(s) = {s-^)logs-s + ^log2n + Oi\s\-') 



and hence 



log- 



1 



,1 



- s) logs + s + 0(1). 



s + 3|a| 



In the strip — | < a < 4, since — | < arg(s + 3|a|) < |, we get 

Re (log ^ ) = ^log s + 3|a| -(a + 3|a|)log 

+ targ(s + 3|a|) + a + 3\a\ + 0(1 
< k2\t\. 

The corollary follows. 



□ 



3. An Asymptotic Formula for T.n{p)<xX'^°'{P) 

Following Kubilius ^ and Landau ^2] we will in this section obtain 
an estimate for J2n(p)<x X^"(P) terms of x and a, by finding zero-free 
regions for the Hecke L-functions L{s, x^") and estimates for the log- 
arithmic derivatives of L{s, x^") in these regions. For the convenience 
of the reader we recall the following two lemmas from Landau's book: 

Lemma 16 (Landau UX,, Satz 374). Let r > 0. Suppose f{s) is 
analytic for \s — so\ < r. Furthermore, suppose 



fis) 



fiso) 



< for \s — sq\ < r 



and 



f{^) 7^ /^'^ k ^ -^ol ^ Re(s) > Re(so 
Then the following holds: 
1) 

-Re 



V/(^o)/ r 

2) // there is a zero p on the line between sq — | and sq ( exclusive), 



then 



V /(-So) / r So- p 
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Lemma 17 (Landau [11 , Satz 225). Let r > 0. Suppose f{s) is 
analytic for |s — So| < r and there satisfies 

Re(/(s)) < M 

Let < p < r. Then for |s — So| < p 



\ns)\ < 



2r 



(r - p)' 



:i\M\ + \f{so)\). 



Remark. From now on, Ci,C2, 
constants (independent of a). 

Theorem 18. In the strip — | < cr < 4 

|L(.,x^")|<Ci(l + |a|)2(l + |t|)^ 



will denote suitably chosen positive 



Proof. From the functional equation we have 

2-2it 

1 



L 



1 



6a 



2n) 



it + 3|a|) 



T{-l + it + 3\a\) 



L[--tt,x"' 



Since 



we get 



6a 



0(1] 



L(--+tt,x' 



< C2 



it + 3|a|) 



r(-i + zt + 3|a| 



Applying twice the functional equation of the Gamma function we get 



i(4 + .«.x-»") 



< C2 



it + 3\a\ 

2 ' ' 



it + 3\a\ 

2 ' ' 



<C3(l+|a|)'(l+|t|)^ 



Furthermore 



0(1). 



6a\ 



m+tt,x 

Consider now the function 

(l + |a|)2(l + s)2 
A(s) is holomorphic in the strip — | < cr < 4, and since 

|L(s,x^")| 



L{s,x 



|A(.)|« 



ii + \a\m + \t\y 



A(s) is bounded on cr = — | and cr = 4 by the above. Furthermore 
it is 0{e'^^) in the whole strip by Corollarv IT31 Thus, by Phragmen- 
Lindelof's Theorem, A(s) is bounded in the whole strip, and the theo- 
rem follows. □ 
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Lemma 19. In the strip 1 < a < 2, 
Proof. 

p 

=(1-3-)-^ n i^-p-T' n i^-p-'T' 

p=l{3) p=2i3) 

= 11(1 n i^-p'T' n i^+p-r' 

P p=l{3) p=2(3) 

where ({s) is the Riemann zeta function, and L{s,x) the Dirichlet L- 
function for the character 

— 1 if n = 2 mod 3 

x{n) = { if n = mod 3 . 
1 if n = 1 mod 3 



Now 



\as)\ < E n-- < 1 + Ju-^du = 1 + < 

n=l { (7-1 (7-1 



L{s,x) = E ^ = 1 u-'d{Y,x{n)} = s ju-^-\Y. X{n))du, 

n=l 



SO smce 



Furthermore, 

X{n 

n=l n<u 2^ n<u 

< 1 for aU a; > 1, 

I^(«,X)I<1 

and the lemma follows. □ 

The next theorem gives zero-free regions for the Hecke L-functions, 
of the same form as can be obtained for the Riemann zeta function and 
the Dirichlet L-functions. 

Theorem 20 (Zero-free Region). There exist positive constants c^,cq 
such that L(s,x^") has no zeros in the region defined by 
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Re 







Zero— free 






region 

/ 




. / 




C6 - 


" \ 


H - 

1 

\ 



Figure 3. A zero-free region for L(s,x^") 



Proof. Logarithmic differentiation of (jH)) yields 



for o" > 1. We recall the definition of 6p on Page ITUl as a unique angle 
in [— |, |). In terms of 6^ we have, for every p and m, 

X*"" log N(f ) _ log N(p) _ ^i^ji 



Ar(p)"^'^ A^(p)' 
Using the inequality 

3 + 4 cos If + cos 2(p = 2(1 + cos ipY > 0, 

we deduce that 

(20) -sm-^^4'4m^]-^4'y^^i^. 



Cxio-) \L{a + it,x^'') J \L{a + 2tt,x 

E (3 + 4 cos (6ma^p - mt log iV(p) 

+ cos (l2ma6'p - 2mtlog A^(p)) ) > 0. 



Let Sq = p + ir, where 1 < p < 2. (p is later to be suitably chosen as 
a function of r.) In the disk \s — So| < | we have, by Theorem 1181 and 
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Lemma ITIll 



6a\ 



6a\ 



(21) 



<c,{l + \a\fil + \t\f L{so,x'n-' 



ci(l + |a|)'(l + |t|)^ 



E 

(a) 



Nia] 



so 



< 



p-1 



(l + |a|)2(l + |r|)^ 



Suppose now that ft + ir is a zero of L{s, x^"), where p — | < < p. 
Then, applying 2) of Lemma with 



M 



C7 



we get 
(22) - Re 



' L'{p + zT, x' 



p-1 



16 



<-log((l + |a|)(l + |r|)) 



1 



- log(p - 1) - 
3 p-p 



Moreover, by 1) of Lemma IT^ 



12a ^1 



(23) 



Re 



' L'{p + 2iT,x 
L{p + 2ir,x^'"') 



<^log((l + |a|)(l + |r|)) 

- ^log(p- 1) +C8. 



As regards ^ we note that, because of the simple pole at 1, 

(24) _ikM 

u(p) p-1 

Now (1201), (E2I), (ESI) and (El imply 

(25) ^<^ + f log((l + |a|)(l + |r|) 
p — 1 3 ^ 



40 



log(p- 1) +Cii. 



p-p p 

We may choose Ce sufficiently large to ensure that for |r| > cg 
401og(1001og2(l + |r|)) + 3cii < 201og2(l + |r|) 
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and thus for all a 7^ 

(26) 401ogfl001og((l + |a|)(l + |r|))) +3cii 



<201og((l + |a|)(l + |r|)). 

Now we put 



1 H -? — i Y for |r| > ce 

1001og(^(l+|a|)(l+|T|)j 

1 H 7 — ^ ^ for |r| < C6 

100 log (l+|a|)(l+|c6|) 



P 



First suppose |r| > cg. Put 

/: = log((l + |a|)(l + |r 

for short. Then ()25|) . multiplied by 3, becomes 
12 

< 980/: + 40 log(100/:) + Scn. 



Thus, by 

< 1000£, 

p-/i 

and hence 

1 12 1 
M < 1 H ^ ^ = 1 ^ 

100/: 1000/: 500/: 

for all eventual zeros /i + it. 
Next suppose |r| < Cg and put 

/:' = log((l + |a|)(l + |c6|)). 

Then 

12 

p -p > 



^ + 80 log ((1 + |a|)(l + |r|)) - 40 log(p - 1) + 3cn 
12 



> 



980/:' + 40 iog(ioo/:') + scn ■ 

But by dSni) (for r = cg) 

401og(100/:') +3cii < 20/:', 

so 

12 

p-p > 



1000/:'' 

and hence 

12 1 12 1 



1000/:' 100/:' looo/:' 500/:' 
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Wc have proven the theorem. 
Lemma 21. On the line (7 = 2, 

logL(s,x''^)<Ci6. 
Proof. U s — 2 + it we have 



27 

□ 



1 



n 



1 - 



N{py 



< 



n (i + m-') 



Thus 



<1 + 1:7V(o)-2 = 1 + Ck(2). 

a 

<\L{s,x''')\<Ck{2), 



1 



i + a(2) 

so the logarithm is bounded. 

Theorem 22. In the region O defined by 



□ 



3> a> i 



1 



ci2log(^(l+|a|)(l+|t|) 
1- 1 



C12log ((l+|a|)(l+|c6|)) 

where Ci2 > 05, u;e /ia?;e 



for \t\ > Cq 

for \t\ < C6 



<Ci3log (1 + |a|)(l + max(|to|,C6 



Proof. Let C14 > C5. For every sq = 2 + ito on the hue Re(s) = 2, let 
be the circle with center at sq and radius 

1 



r — < 



1 + 



1 + 



C14 log((l+|a|)(l+|to|) 
1 



C14 log (^(l+|a|)(l+|c6|)) 

We may freely assume that ce > |, so that 

(1 + |a|)(l + C6) > 3 
Then we have for s = u + ^^ in CgQ : 



for 1 > Cg 

for \to\ < ce 



log((l + |a|)(l+|t|)j <log((l + |a|)(3+|to|)^ 

<log ((l + |a|)(l+|io|))+log3 

Thus we get 

'21ogf(l + |a|)(l + |io|)) if|io|>C6 
21og((l + |a|)(l + |c6|)) if|to|<C6 



(27) log((l + |a|)(l + 



< 
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Now we want to use Lemma fTTI with f{s) = logL{s,x^°'), which is 
analytic in the zero-free region of Theorem 1201 and thus analytic in 
C 

Also, by Theorem ITHl 

\L{s,x'n\<ci{l + \a\fil + \t\f 

in Cso, so by (j2I|) 

Re(/(s)) = Re(logL(s,x''^)) 

<41og((l + |a|)(l + |t|)) +logci 

< 4 log (^(1 + |a|)(l + max(|to|,C6))) + C15. 

Set 

1 H J — i ^ for |to| > C6 

C12 log (l+|a|){l+|to|), 



1 + ^ T for |to| < C6 ' 

C12 log(^(l+|a|){l+|c6|)j 



where C12 > C14. Now Lemma IT7| with M = 41og((l + |ci|)(l + 
max(|to|, Cq))) + Ci5, and Lemma 1^ imply that in the disk |s — sqI < p 



we have 



L'{s, x'"") ^ _2r_^ |^(^ ^ 1^1^^^ ^ max(|to|, ce))) + c,, + |/(.o)|) 

41og((l + |a|)(l + max(|to|,C6))j +016). 



(r — pY 
2r 



< 



(r — pY 
But 



F^pF " f 1 _ 1 f l^^ + (1 + max(|to|, C6 

\Cl4 C\2 ) 

- ( 1 ^ 1 f (^^ + '""'Kl + max(|to|, ce))) . 

VCl4 C12 J 



The statement of the theorem follows. □ 
Definition. Let 

This series is easily seen to be absolutely convergent for a > 1, and 
uniformly convergent for a > 1 + 6, so K{s, x^") is analytic for a > 1. 
We even have: 
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Lemma 23. K{s,x^"') is analytic in the region fl of Theorem\EE and 
satisfies 

' Ci5log=^f(l + |a|)(l + |t|)) >C6 
Ci5log=^((l + |a|)(l + |c6|)) Zf\t\ <C6 



\K{s,xn\< 

there. 

Proof. By Proposition EJ we have, for cr > 1 

But the logarithmic derivative is analytic in VL, and the sum on the 
right is absolutely convergent for ex > ^ and uniformly convergent for 
cr > I + 5, (5 > 0, since 

logiV(p) _^ logiV(p) 



p n^=2 iV(p)™(l+^) V iV(p)i+^(iV(p)l+^ - 1) 

logp^ 



<2E- 



p 



oo 



- 1) 



^ ^ log 
< 4 > ^— — ^- < oo. 

Thus, constitutes an analytic continuation of K{s, x^"") to VL (since 
clearly VL lies to the right of the line a = | ) . From (|^Hj) it is also clear 
that, in f2. 



< Cl7. 



□ 



L(s,x6'^) 

Thus in f2 we have 

<Ci3W ((l + |a|)(l+max(|to|,C6))) +ci7 

< Cislog^ {{I + |a|)(l + max(|to|,C6) 
Lemma 24. 

2+ioo 

/■ —ds -1° «/0<x<l 
J s'^ I 27ri log a; if x > 1 

2-ioo ^ 

Proof. Let x > be fixed. The function ^ is analytic in the whole 
plane, except for a double pole in the point s = with residue logx. 
Assume first that < x < 1. Using the integration contour of Figure 
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2 + iR 




13] we have by Cauchy's Theorem 

2+iR 

X 



X 



—ds + / —ds = 0. 



2-iR 



But 



x^ 



ds 



rtR 



< vr 



X 



SO letting i? — ^ oo we get 



2+ioo 



X" 



ds = 0. 



2— ioo 



Assume instead that x > 1. We then use the contour of Figured Since 
(if we take R> 2) the pole lies inside the contour we get by Cauchy's 
Theorem 

fx". /■ X% ^ . , 
j —ds + J — ds = 2m log X. 



2-iR 



But 



x^ 



ds 



< TT 



X 



SO again we let i? ^ oo, yielding 

2+ioo 



fx, , 
J —ds = 2m logx. 



2—ioo 



□ 
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Now we are ready to estimate sums of characters. We pave the 
way to Y.x{p) by estimating two weighted sums, with the weights 
log A^(p) log and logA^(p), respectively. For the first one, we will 
use the integral of Lemma 1^ to pick out the portion of the L-series for 
which N{p) < X. We will then shift the integration path a bit to the 
left, close to the zero- free region. 



Theorem 25. For x > 1 



X^"(p)logiV(p)log— ^ < ^''^°g(l+l''^)+^^^log^(l + |a| 

N{p)<x ^^^) 

Proof. By Lemma 1^ 

2+joo 2+joo / X 



(29) 2-ioo P 2-ioo 

= y: x''^(p)iogiv(p)iog-^, 

which is the sum we want to approximate. Now let u be the curve 
defined by 

1 7 — Y for |t| > ce 

ci2log{(l+|a|){l+|t|) 

1 ^ ' for|t|<C6 ' -^<t<^- 

ci2 log(^(l+|a|)(l+|c6|)j 

I claim that 

2+ioo ^ ^ 

(30) J ^K{s, x"')ds = J ^K{s, x"')ds. 

2-ioo <^ 

To see this, consider for large T the contour defined in the following 
manner (see Figure E)): 

from 2 — iT to 2 + zT in a straight line, 

from 2 + zT to 1 ^ + zT in a straight line, 

C12L ' 

from 1 + iT to 1 ^ — iT along u 

C12C C12C ° 

and from 1 ^ — to 2 — zT in a straight line. 



(Here we have written for short C = log y{l + \a\){l + |T|)j.) 
Since the integrand is analytic inside and on ■jt, by Cauchy's Theorem 
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T 
1 - 




- — 










It 






/ 




C6- 


I 






-C6- 




H 1 

1 

\ 




-T- 









Figure 6. 



= J^K{s,x'nds 



IT 



(31) 



/ 2+iT 



NT 1- 



2-iT 2+iT i_ 



-iT 



1- 



2-iT 



\ 



-iT 



But by Lemma 



2±iT 



1 ^±iT 



<|^Ci5log3((l + |a|)(l + T) 



so letting T — > oo in (p?T|) the horisontal integrals vanish, implying 
Thus we need to approximate the integral along u. Now, for an 
arbitrary r > cg, we have by Lemma 1221 



C6 1-- 



J ^K{s,x'nds^ J log='((l + |a|)(l + C6))cit 



T 00\ \ 



X ci2log{{l+|a|)(l + t)) 



log^ ((1 + |a|)(l + t))c/t 



\C6 T 
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The right side is 

< 2;l-ci2log{{l+|a|)(l + C6)) log^(l + |a|) 

1 I 7^og'{{l + \a\){l+t)) 

-\- X ci2log((l + |a|){l + T)) / V L(]^f 

J 



T 

Putting together the first two terms we get 

log ^ n 

< xe <=i2i°g((i+ki)(i+-)) log (1 + \a\) 
+ -logVlog^(l + \a\). 

T 

Putting r = eV^^, we get 

But for large x, 

^ log log X - Jhgx < -C2o\/logx < -C20: 



2 ^ ^ log(l + \a\) + y/\ogx 

so, putting C18 = min(ci9, C20) and recalling and pUj). we deduce 



that 



X^" (p) log N{p) log <^ xe ''''i°g(i+i-i)+v^ log^(l + |a|). 

Nip)<x ^^^) 

□ 

Theorem 26. For x > 1 

log X 



J2 X^^ip) log A^(p) -C xe iog(i+ki)+0^ log3(l + \a\). 

N{p)<x 

Proof. Set for short 



S = 6{x) = e ' iog(i+ia|)+v^. 
With X replaced by (1 + 6)x, Theorem gives 

E X^"(P)logiV(p)log^i±^ 

N{p)<{l+S)x ^^IPJ 

(32) „ log(l + f):r 

<(l + (5)xe iog(i+ki)+Vi°g{i+^)- log^(l + |a|) 
< 5^a;log^(l + |a|). 
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We will now split the sum on the left side in two parts. First, again 
using Theorem 123 we have 

(33) 

y: x'"(p)iogiv(p)iog^^^±^=iog(i+5) y: x'"(p) log iv(p) 

N(p)<x N(p)<x 

N{p)<x ^^IPJ 

= log(l + (5) Y x'"(p)logiV(p) 

N(p)<x 

+ 0{6^x\og'il + \a\)). 

Secondly, 
(34) 

Y X'"(P) logiV(p) log ^ii-^ « 5xlog((l + 6)x) log(l + 6) 

x<N{p)<{l+5)x ^^IPJ 

■C S'^x logx, 

since the number of terms in this sum is 0{6x). By (|33|) we have 
log(l + 5) Y x'^iP)^ogN{p)= Y x''^(p)logiV(p)log^y^ 



N{p)<x N{p)<x 

+ 0{S^ x\og^{l 



Y X''^(P) log iV(p) log + 



Nip) 

Y x'"(p)logiV(p)log 



N{p)<{l+5)x 

il + 6)x 



x<N{p)<{l+S)x 

+ 0(S^x\og\l + \a\) 



N{p) 



By dSl and dni, this is 



0(5^xlog^(l + \a\)j+ 0{6\ logx) + 0[6^xlog%l + \a 
0((5^xlogxlog^(l + \a\) 
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whence 



^2 

J2 x'"(p)logiV(p) < xlogxlog=^(l + |a|) 

Trrt.. log(l + 0) 



N{p)<x 



<^ 6x logo; log (1 + \a\) 



'■^'^18 ]ogJ£ +log log X , n, , ,. 

a;e 2 «iog(i+|a|)+vi^ log^(l + |a|) 

log X 

3;g ''^Sog{l+|a|) + v^ l0g^(l + |a|). 



□ 



The final result of this section now follows easily: 
Theorem 27. For x > I 



log X 

X^^ip) =< xe~ log^(l + \a\). 



N{p)<x 



Proof. Let 



N(p)<x 



1 \ ^ i}{x) 



We employ partial summation: 

N{p)<x 2<m<x lOgm 

By Theorem OBI this is 



log m /I 1 

< E '^^i°g(i+i''i)+ViS^ log^(l + |a|) 

2<m<x 



m log(m + 1)/ log([x] + 1^ 



logm log(m + 1"; 
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-C21 '-^ 



But since xe iog{i+\a\)+^n^ jg monotone increasing for sufficiently 
large x, 

'1 1 \ 



loK m 

-C21- 



^ me iog{i+|a|)+viogm log^(l + \a\) 



2<m<x 



\ogm log(m + 1) 



< a;e~''^''°*5(i+i'^!)+Vi°g- log^(l + 



2<m<x\l0g"^ log(m+l) 



log X /I X 

< xe log^(l + |a|) 



^log2 log([x] + l) 



■C xe '"sCi+i-D+v^ log^(l + |a|). 



The theorem follows. □ 

4. Equidistribution of lattice points on circles 
4.1. A Lemma on multiplicative functions. 
Definition. A function / : N ^ C is called multiplicative if 

f(mn) = f{m)f{n) 

whenever {m,n) = 1. 

We will prove a form of the so called Halberstam-Richert inequality, 
giving a bound for Y.n<x /('^) via a bound for Y.p<x This version 
appears in Katai [19j. 

Lemma 28. Let f{n) be a nonnegative multiplicative function, and 
assume that 

fip'^) = 0{a) 
for every prime power p"", a > 1. Then we have 

Proof. Let 

Aix) := E f{n), B{x) := ^ /H ^ogn. 

n<x n<x 

Then we have 

i3(x) = ^/H^logg° 

n<x qO'Wn 

= E /(/^)/(01ogg", 

q^hKx 
(<?-,h) = l 



HEXAGONAL LATTICE POINTS ON CIRCLES 37 

since f is multiplicative. Hence 

(35) B{x)=J2f{h) Y: /(g")logg^ 

h '3"<f 

But for an arbitrary y > 2 we get, using the assumption on /, 
Y: /(g")logg"« ^ a2logg = ^a2 ^ logg. 

By the Prime Number Theorem (see for example §18) 

Y logg = 0(y«), 

X 

q<ya 

and thus 

Y fm log g° « E =y + T. ^'y~- ■ 

q"<y a a>l 

For every a in the sum we have 2" < y, so 

E"'2/"< E a^y^ < E (r^) 2/^ < 2/^ log' I/, 

since the number of terms in the sum is 0{logy). Thus 

E /(?°) log 9" = y + 0{y-^ Wy) = 0{y). 

q°'<y 

Inserting this into ()35|) yields 

(36) 5(x)«E/W^ = ^E^- 

h<.x h<x 

But since / is nonnegative and multiplicative, we have 

h<x p<x \ y y / 



p< 

We claim that 



p<x\ P J \ P P J 



< oo. 



Indeed, this product converges absolutely if and only if the series 

p\p p^ J 
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converges absolutely, and by the assumption on / this series is 

P k=2 f 

Thus 

so, recalling (jHUj) . we get 

(37) B{x)<^xexplY.^ 

Now we can make a crude estimate of A(x): by partial Stieltjes inte- 
gration we have (assume without loss of generality that x > 2) 
(38) 

A(x) = 1+ y fin) = 1+1 r^dBit) = 1 + ^^^+ / -^^B(t)dt. 



2<n<x 



2- 



But by Merten's Theorem (see for example Ch. I.l, Theorem 9), 
(|H7j) implies 



B{x) <^ a;exp ■{ E ~ }■ ^ a;(loga;)'^ 



and hence 



X X 

f ^^B{t)dt<^ [{\ogtr-'dt<^ 

J tloK t J 



X if c < 2, 

a;(logt)'=~i ifc>2. 

Putting these estimates into yields 
(39) A{x) < 

for an arbitrary e > 0. Now 

A{x)-A{V^)= y fi^)<A- E /Wlogn 

lOET X 

•/x<n<x ° s/x<n<x 

<-^5(x)«-^exp|5:M 
logx logx p 

But by (jni, 

^(^) = 0(x^/2+.)^ 
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SO this term is negligible, and the lemma follows. □ 

4.2. Proof of Theorem [31 We recall the definition of the exponential 
sums S{n,A) on Page IHl and define functions fA^n) in the following 

way: 

Definition. Let 



fA{n) 



\S{n,A)\ 
6 



These will prove to be multiplicative functions of n, thus allowing us 
to make use of Lemma I 



Lemma 29. fA{n) is a multiplicative function. 

Proof. Suppose that n has the prime factorization of (jH)) on page |H1 If 
we exclude the trivial case where one or more of the /5j are odd, and 
consequently /^(gf*) = = /^(n), we get 

m=Oji=l ifc=l 

Since 

5 



l(a7r/2+(ai-2ii)epj+...+{Qfc-2ifc)epj^) 



(40) 

we see that 
fA{n) = 



m=0 



6 if A = (mod 6), 
otherwise, 



\Sin,A)\ 
6 



lA 



E---Ee 

ii=i ife=i 

^ giA(ai-2ji)6»p^ 
ii=l 

fAiPl)---fA{Pk), 



((ai-2ji) 



K-2i 



ife=i 



so /a is multiplicative. 



□ 



Furthermore, by (|40|) . /a = when A ^ (mod 6), and in this case 
all results are trivial, so we henceforth assume that 6 | A and substitute 
A with 6a, a G Z. 

We examine the values of f^a for prime powers: 

. /6a(3") = 1. 

• If g = 2 (mod 3) we have 



/6 



1 if a is even, 
if a is odd. 
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If p = 1 (mod 3) we have 



j=0 



^6ia(a—2j)9p 



<a + l. 



Thus fsa clearly satisfies the hypothesis of Lemma |2H1 Moreover, we 
have in particular 



(41) 



2|cos6a^^p| if p = 1 (mod 3), 
/6a(p) = <|0 ifp = 2 (mod 3), 

1 if p = 3. 



Now we need to calculate Y,p<x • From Fourier analysis we 



recall Fejer's Theorem (^Hl! Th. 1.5), stating that if is a continuous 
function on [0,27r], then the arithmetic means of the partial sums of 
the Fourier series of g converge to g uniformly on [0, 27r]. Thus for each 
e > there exist k and ao, . . . ,ak such that for every x G [0, 27r] 



cosx 



am COS mx 

m=0 



< 6. 



There are only cosine-terms in the Fourier series, since the function 
I cosxl is even. Moreover 



2n 



ao 



— / I cosxMx = — . 



27r 



TT 



By (HH) we get 



E 



p 



+ E 



2 1 cos6a6'„ 



p<x 
p=l (3) 



P 



L 2am cos QamOr. 



(42) 



p<x 
P=l (3) 

•2 



v,m=0 



P 



2e_ 
P , 



IT 



E - 

p<x P 
p=l (3) 



+ E '^'i 



y-v 2 COS QamOp 

p<x P 
\p=l (3) 



In order to estimate the sum on the right we will reformulate Theorem 
123 a bit: 
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Proposition 30. If x > 1, a ^ and \am\ = 0(ev^^) we have 

p<x 
p=l{3) 

Proof. We take a closer look at the sum 

N{p)<x 

If p is a prime ideal, then either A^(p) = p, where p is a split prime, or 
A^(p) = q"^, where q is an inert prime, or A^(p) = 3. In the first case we 
have (p) = pp = (7rp)(7rp). Thus, if p = 1 (mod 3) we have 

E ^'"'"(P) = X'"'"(vrp) + X'""^(vrp) = 2cos6am^,. 

N{p)=p 

If g = 2 (mod 3) we have 

iV(p)=q2 

Finally 

E X'""^(P) = X'""^(vr3) = (-I)"'". 

JV(p)=3 

We conclude that 

(43) E X^'""(P) = E 2cos6am^p+ E 1 ± 1- 

N{p)<x p<x q'i<x 

P=l(3) 9=2(3) 

Moreover, if 1 + \am\ < Ce"\/'°^^ we have 
logx 



-= > Cyiogx, < C" < 1, 
log(l + \am\) + vlogx * 



so by Theorem 1771 



log X 

-C21- 



N{p)<x 

3 



<a;e-"^Vi°g^(loga;)5 
Thus by (021) 

E 2 cos6am^p < xe'^^^^v^ + < xe'''^^^^, 

p<x 

p=m 

and the proposition follows. □ 
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Note the restriction on a. A simple partial integration now yields 
the estimate we want: 

Proposition 31. If x > 1, a and \am\ = 0(eV^^°^) we have 
(44) J2 1^2^^^ = 0(1). 

p<x P 
P=l{3) 

Proof. By Proposition IHUl 

p<x P 5_ I- P<t ) 

p=l{3) p=l{3) 

X 

< 0(1) + e-'^^^v^ + j t-^e~'^^V^'dt. 

5 

Now, by the change of variables u = y/\ogt, we have 

< / t-^e-^2*v^dt < / r^e-'^^^^'dt = f 2ue"^2*"dn = — , 

^24 



5 

and thus 



^ 2cos6amgp « e'^^^v^ + 0(1) = 0(1). 

p<x P 

p=m 



□ 



From the Prime Number Theorem for Arithmetic Progressions ( ^3] , 
Ch. 20) it easily follows that 

(45) J2 - = ^loglogx + 0(l). 

t^x P 2 

p=m 

Inserting (jSj) and into (jl21), we get 

fM / 1 , ^2 

p<x P 



51 P 3 Vtt / V i<m<fc' 7 



2 

+ e ) log log X + Oe(l), 



vr 



uniformly in a, provided a 7^ and a = Oe(ev'°^). 
Now, using Lemma EHl we conclude that 

y)/6a(^) < ^j-^ exp ( f- + £:) loglogx + Oe(l)| <x(logx)- 
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which proves Theorem El □ 

5. The Distribution of Hexagonal Primes in Sectors. Bad 

Circles. 

Aheady in 1903 Landau J13J proved the so called Prime Ideal The- 
orem, giving an asymptotic formula for the number of prime ideals 
with norm < x in an arbitrary number field, then with the error term 
0{xe~^^°^^'*^^^^). The version we give here comes from First we 

define the logarithmic integral: 

Definition. For x > 2 we define 




Theorem 32 (Prime Ideal Theorem). For the number of prime ideals 
in a number field K with norm < x, we have 

7r^(a;) = ^ 1 = Li(x) + 0(xe~^v^), 

N[p)<x 

where n is the degree of K , and h is a positive constant, independent 
ofK. 

We will refine Theorem in the particular case of the hexagonal 
number field K = Q(v^— 3) to a result measuring the number of prime 
ideals p < D such that A^(p) < x and such that 9p lies in a specified 
interval y92]- Again we follow Kubilius [H]. 

To this end we will apply the results on the previous section to a 
Fourier series. We will use a lemma of Vinogradov to prove that the 
characteristic function of the interval [ipi, 992] can be approximated by 
functions with well-behaved Fourier coefficients. 

Lemma 33 ([ZI], Ch.l, Lemma 12). Let r be a positive integer, and 
let a,P,A be real numbers satisfying 

0<A<1, A</3-a<l-A. 

Then there exists a periodic function ipi^), with period 1, satisfying 

(1) ip{x) = 1 in the interval a + ^<x</3 — ^, 

(2) ^jj{x) = in the interval /3 + ^<x<l + a — ^, 

(3) < ^p{x) < 1 in the remainder of the interval a — ^ < x < 
1 + a - f , 
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(4) ipi^x) has an expansion in Fourier series of the form 

oo 

ip{x) — (3 — a + ^ {am cos 27rmx + bm sin 2Trmx), 

m=l 

where 



1 (^m 1 ) 


\bm\ 


l^ml ) 


\bm\ 



Remark. If we instead express the function tp{x) above in a Fourier 
series of the form 



27rmx 



then we have 



CQ^(3-a, Cn^^{an-ibn), c_n = ^(a„ + i6„), (n > 1) 



so that 



|c„| < 2(7r|n|)-S 
2 / r 



c„. < 



7r\n\ \7r|n|A 



Lemma 34. Let S > and suppose 25 < </72 — V'l < f — 25. Then there 
exist ^-periodic functions f{(p) and f_{(p) such that 

(1) 7(<^) = 1 «/<^l < <^ < ¥'2 

< /('P) < ^ in the rest of the interval (fii — S < (fi < 

(2) = l «/</^i + 5<¥^<¥^2-5 

/(99) = 0z/^2<V5<f + V'i 

< < 1 in the rest of the interval (fi<(p<^ + (fii. 



(3) ^/ 



n=— oo n=— oo 

i/ien we have 

-3, _ 1 _ 1 

Oo = -(</^2 - + 5), a„«:--, a„ < — 

TT \n\ d\n\ 

3, 1 1 



2 



TT |n| d\n\ 



2 ' 
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ip2 



Figure 7. 



Proof. This follows immediately from Lemma IH^ if we take x = ^ip and 
r = 1, setting for /, 



and for / 



3 3~„3 
a = -(pi - —5, /3 = -ip2 + — 5, A = -5, 



3 3_„3 
a = -(pi + —6, (3 = -(p2 - — 5, A = -S. 

IT ZTT TT /TT TT 



Definition. Let — ^ < v^i < v^2 < f - Then we define 



iV(p)<x 

i/ji<ep<</32 



□ 



We now prove the theorem about the distribution of prime ideals in 
circle sectors. 

Theorem 35. We have 



TT 



[^,,^,](a;) = -(^2 - ^i)Li(x) + 0(xe-^-V'°s^). 



TC 



Proof. Define the functions /(v^) and fi^p) as in Lemma IMl with 
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Then we get 



oo 



N{p)<x N{p)<x N{p)<xn=-co 

(pi<9p<ip2 

= ao7r(x) + ^a„ Yl x'"(p) 

n^O \N{p)<x 



Thus, by Theorems ESI and |2 

3 / 



C27ylOgZ^ 



(46) 



-C21- 



+ \ Y1 l«n|xe iog{i+|n|)+v^ log (1 + |n|) 



Analogously we deduce 



vr[^„^.](x)> E m)= E E ^nX'"(p) 

Af(p)<a:: Ar{p)<a: ™=-oo 

3 



(47) =-{^2~^i-S)(Li{x)+0{xe- 



C27\/ lOgX^ 



-C21- ^ 



We examine the sum on the right side of ()46|) . By the bounds on |a„| 
in Lemma OH we get, if we split the sum in two parts: 



E 



xe 



-C21 



1 og { 1 + I n I ) + y/log~r 



log^(l + \n\) 



^ log^n -C21 , ^ log^n 



l<n<5-2 



For the first part we note that 



log(l + 5 ^) < log5 ^ < V^logx, 
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and thus 

E 

l<n<<5-2 



\og^n -C21 



log X 



n 



3/'A-2\ „-C28 Vl°g^ 



E 



l<n<(5-2 



n 



< e-'=^«vi^log^(5-2) 
= e-'=^«v^(2c26V^)' 



For the second part we have 



n>(5- 



Obviously we have the exact same bounds for the corresponding sum 
in (jlTJ containing a„. Thus (gHl) and (gH) yield 



7r[^i,V2](a;) - -{^2 - ¥'i)Li(x) 



•C xe 



-C25 V logX 



and we are done. 



□ 



Since there is a one-to-one correspondence between prime ideals and 
hexagonal primes (that is, primes of the number ring O) in the angular 
interval [— |, |), we get 

Corollary 36. The number of hexagonal primes in the circle sector 
{z; \z\ < Vx, < argz < ip2}, 



TC IT 

where - - < < ^2 < - 
b b 



IS 



-{ip2 - (/^i)Li(x) + 0(xe-^25Vi°g-). 



We also have 



Corollary 37. The same estimate holds if we consider only non-real 
primes. 

Proof. The number of real hexagonal primes in the sector described 
above is obviously 

□ 



Theorem E] also follows directly from Theorem |HS1 
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Proof of Theorem^ By Theorem 1221 and Corollary IHEl we have, for 

-i < V'l < V2 < f 

7r(x) (l + o(l))Li(x) vr^^ 

This is what is required in the definition of equidistribution on Page 

IS □ 

We will now show that it is possible to construct arbitrarily "bad" 
circles. We follow Cilleruelo |22j . 

Theorem 38. For every e > and every A; G N there exists n G N 
such that the circle with radius \fn centered at the origin has more than 
k lattice points, all of which are concentrated on the six arcs 

{V7^e*('^i+'^); |(^| < 4, u = 0,l,...,5. 

Proof. Let e and k be fixed. Choose an integer m such that 

log k — log 6 

m > — , 

log 2 

and let < 5 < ^. Then by Corollary |3] we can find m different primes 
Pi,...,Pm such that 6 < 6p^ < ^. Set 

n = pi---p^. 

Then the solutions a G D to the equation n = aa all have the form 



In each case 



Q, ^ ^^i{±ep^±...±ep^)+iuf 



±6p^ ± . . . ±6p^ \ < e, 



and by ((Zj) the number of solutions is 

rQ{n)=6-2"'>k, 
which proves the result. □ 

An example clarifies the method: 
Example. Put 

n = 7983607 = 157 -211 -241. 

We have 

157 = 12^ + 12 ■ 1 + 1^, ^157 ^ 0.0692 < — , 



36' 

211 = 142 + 14-1 + 12, ^211 ~ 0.0597 < —, 

36 

241 = 15^ + 15 ■ 1 + 1^, ^241 ~ 0.0558 < — , 

36 
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SO that 

TT 

Thus the lattice points on the circle with radius ^/n, centered at the 
origin, all have arguments lying, modulo |, between and ^. More- 
over, this circle has 6 ■ 2^ = 48 points. (See Figure |H1) 

6. A Further Measure of Equidistribution 

We will in this section give a result that perhaps better justifies the 
statement of equidistribution on average. In the case of true equidis- 
tribution, the ratio of points lying in a specific subinterval would be 
approximately equal to the ratio between the length of the subinterval 
and the length of the whole interval. The discrepancy A{n) measures 
how far off this approximation is. 

Definition. Let 




-3000 -2000 -1000 1000 2000 3000 



Figure 8. Lattice points on the circle with radius r = V7983607. 
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We will show that the discrepancy is small for all but a few values 
of n, as n grows large. Obviously we must specify the meaning of "a 
few" - from Section ^ we know that in some sense only "a few" circles 
with radius -Jn, n G N, have any lattice points at all (this happens if 
and only if all prime factors q = 2 (mod 3) occur in even powers in n). 

We recall that Q is the quadratic form defined by 

Q{x,y) = + xy + y'^, 

and that rQ{n) is the number of representations of the integer n by the 
form Q. We put 

TZq{x) := {n<x ; rgin) ^ 0}, 

BQix) := |7^Q(x)|. 
The asymptotic formula for Bq{x) was found by Landau 

ts a consta 

Bq{x) ~ b- 

V^ogx 

From now on let r(n) = rQ{n). Our result will take the form: 

Theorem 40. For almost alln G TZq{x), that is, with the exception of 
o{Bq{x)) of them, we have 

A(?7.) < rqin)"^, 

In proving this we follow Katai and Kornyei , who gave the anal- 
ogous result for the square lattice, as did Erdos and Hall ilOJ. We will 
need a result of Erdos and Turan [23j concerning the discrepancy: 

Lemma 41. Let Lpi, . . . , (p^^ G M. Put 

1 ^ 
Zu = — Y^ e''"^^ 



Theorem 39. There exists a constant b > such that 

X 



Then for an arbitrary T > we have 



sup 

0<a</3<27r 



1 l_ ~^ 

AT 2^ 



a<(pj<fi 



T 



1 V- 



fc=i 

Proof. See [23!, Th. III. □ 

Proof of Theorem For < 7 < 1, put 

C,{x) = E A{n)r{nr, 



n<x 
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and 

M^(x) = #{n < X ; A{n) > r{ny^}. 

Our aim will be to prove that C^{x) = o{Bq{x)), since then it follows 
that My(x) = o{Bq{x)). By Lemma 1^ we have for an arbitrary T > 

1 T 1 |^(^^)| 



and thus 



T .^r; A 



n<x n<.x A=l \ ' 

n<x A=l n<x 

«^ + ETE6^^A(n), 

A=l ^ n<x 



where 



gA{n) = Q-'\S{n,A)\r{ny-^ = fA{n) 



' r{n) 
~6~ 



7-1 



gA{.n) is easily seen to be a multiplicative function of n, and since for 
primes p we have 

gA{p) = 2^-Va(p), 9A{pn <{a + ly = 0(a), 
another application of Lemma OHI on Page EHl yields 

E 9A{n) < x(\ogx)-^'-\ 

n<x 

for an arbitrary e > 0. Putting T = [logs] + 1 now yields 

C^{x) <^ h a;(logx)~^^~Mogloga;. 

logx 

:hoose 7 < ^ - 
' log 2 

we get 



If we choose 7 < — 1 and e sufficiently small, so that ~~ + ^ < |) 



C,{x)=ol-^)=oiBQix)), 
\^\ogxJ 

which completes the proof. □ 
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